Hydrodynamics of self-propelled hard rods 



o 
o 

(N 
< 



O 



O 

o 



o 

00 

O 
l> 
O 



X 



Aparna Baskaran and M. Cristina Marchctti 
Physics Department, Syracuse University, Syracuse NY 13244 
(Dated: February 1, 2008) 

Motivated by recent simulations and by experiments on aggregation of gliding bacteria, we study 
a model of the collective dynamics of self-propelled hard rods on a substrate in two dimensions. 
The rods have finite size, interact via excluded volume and their dynamics is overdamped by the 
interaction with the substrate. Starting from a microscopic model with non-thermal noise sources, 
a continuum description of the system is derived. The hydrodynamic equations are then used to 
characterize the possible steady states of the systems and their stability as a function of the particles 
packing fraction and the speed of self propulsion. 

PACS numbers: 87.18.Ed, 47.54.-r, 05.65.-|-b 



I. INTRODUCTION 

Self propelled particles draw energy from internal or 
external sources and dissipate this energy by moving 
through the medium they inhabit. Such a definition en- 
compasses a wide class of systems such as fish schools, 
bacterial colonics and monolayers of vibrated granular 
rods. In all of these the energy input that maintains the 
system out of equilibrium is on each unit, rather than at 
the boundaries as in more conventional non-equilibrium 
situations . A striking phenomenon exhibited by these 
systems is flocking, the emergence of a coherently mov- 
ing body of self propelled entities of size much larger than 
the length scale of the intcr-particlc interaction. 

Extensive theoretical effort has been devoted to un- 
derstanding these non-equilibrium systems. Two distinct 
approaches have been used. First, starting with the sem- 
inal work by Vicsek a number of numerical studies 
of simplified rule-based model systems have been carried 
out. A second approach has been to write down generic 
continuum theories based on symmetry considerations 
0. These methods are well developed for equilibrium 
systems, where general properties such as the fluctuation- 
dissipation theorem yield strong constraints on the pa- 
rameters of the dynamical model [1, [^. They are also 
very useful for non-equilibrium systems, although in this 
case the parameters in the equations remain largely un- 
determined. Both approaches have exposed the dramati- 
cally different nature of order-disorder transitions and of 
the fluctuations in the various phases exhibited by self- 
propelled systems when compared to their equilibrium 
counterparts. Collections of self-propelled units can ex- 
hibit long range order in two dimensions Q , in sharp con- 
trast to equilibrium systems, where the continuous rota- 
tional symmetry cannot spontaneously be broken in sys- 
tems with short-range interactions @ . Recent numerical 
work has shown that the Viscek model and other closely 
related models with polar aligning interactions exhibit a 
first order (discontinuous) transition from a disordered 
to an ordered state in two dimensions when the density 
and the noise amplitude in the system are varied • On 
the other hand, when this model is modified to induce 
strictly nematic order in the system, the phase transi- 



tion becomes continuous [sj]. In addition, large number 
fluctuations have been predicted and observed in the ho- 
mogeneous states of these driven systems [1, [l^, [Hi • 

Non-equilibrium statistical mechanics is a powerful 
tool that can be brought to bear to describe this rich class 
of systems. One class of questions pertains to the iden- 
tification of the underlying collective mechanisms that 
give rise to emergent behavior. Comparison of theoreti- 
cal and numerical studies of model systems can provide 
insight into the physical origin of the various phenomena 
captured by the generic continuum theories. The second 
class of questions is associated with identifying minimal 
microscopic models capable of accounting for a given ob- 
served phenomenon. This is important when trying to 
generalize concepts developed in the context of animal 
group behavior to artificial systems considered in collec- 
tive robotics [iJl ■ An exarnple of this approach is the re- 
cent work by Bertin et al. [l3| who used non-equilibrium 
statistical mechanics to derive continuum hydrodynamic 
equations for the polar Viscek model, with parameters 
given explicitly in terms of parameters of the microscopic 
model. 

In this work, we study the collective behavior of self 
propelled hard rods with excluded volume interactions, 
moving on a substrate in two dimensions. Part of our mo- 
tivation comes from recent simulations of this systems by 
Peruani ct al. [ij]. Unlike earlier work [H, [11], our the- 
ory accounts for the extended shape of the self-propelled 
units. Wc consider a minimal model, with the aim of 
identifying the simplest physical mechanisms that can 
lead to self-organization in self-propelled systems with 
short-range interactions. The non-equilibrium statistical 
mechanics for this model is developed systematically and 
results in coarse-grained hydrodynamic equations for the 
slow variables of the system. The macroscopic equations 
are then used to characterize the homogeneous steady 
states and their stability as a function of particle density 
and self propulsion velocity. This analysis examines the 
effect of flow from self propulsion on the intrinsic entropic 
order that arises from the finite size of the rods. Our work 
provides a microscopic basis for some of the phenomena 
predicted in active nematic via generic continuum theo- 
ries [l^ and observed in simulations d, [3, [2^ . 
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A real system that may be described by our minimal 
model is myxobacteria on a substrate under starvation 
conditions. These bacteria are rod-like in shape, with 
an aspect ratio of about 7, and retain their shape under 
movement. Under starvation conditions the C-signalling 
mechanism in these bacteria that give rise to reversal of 
the direction of motion is suppressed and alignment is 
thought to be brought about by steric interactions [l6j . 

The layout of the paper is as follows. First, the mi- 
croscopic model is introduced and the non equilibrium 
statistical mechanics developed. In a low density approx- 
imation, this is described by a Smoluchowski equation. 
From this starting point, the hydrodynamic equations as- 
sociated with the slow variables of the system, namely the 
density, polarization (related to the self propulsion flow 
field) and ncmatic order parameter, are derived. Next, 
the homogeneous steady states are identified and their 
linear stability with respect to spatial fluctuations char- 
acterized. This analysis is used to interpret the physical 
mechanism behind the emergence of inhomogeneities in 
the nematic state observed in the numerical work of Ref. 
[l^ and the large number fluctuations predicted in . 



II. MICROSCOPIC MODEL 

As a minimal physical model for a system of self pro- 
pelled particles, we consider a collection of rigid rods 
of length i >> b, with b the rod diameter, moving on 
a substrate in two dimensions (2d). Each rod is "self- 
propelled" in that it has a force F acting on its center of 
mass and oriented along its long axis. The rods interact 
only via excluded volume interaction and their dynam- 
ics is overdamped. The exchange of momentum with the 
substrate due to friction provides the physical mechanism 
that results in all the rods on average having a constant 
velocity of magnitude vq along their long axis. At high 
packing fractions, the excluded volume interaction gives 
rise to orientational order, as in passive nematic liquid 
crystals. 

Note that a flock of self propelled particles is charac- 
terized by three properties : 1) centering - the ability of 
a collection of such particles to stay together, 2) velocity 
matching - all the particles are on average characterized 
by the same velocity, and 3) alignment - the particles 
are asymmetric and exhibit ordered states with either 
nematic or polar symmetry p7| . In our model, the rods 
have no attractive interactions and hence will not form a 
cohesive flock. We assume that the system is at a fixed 
homogeneous density po, which is maintained by confine- 
ment at the boundary of the system. Alignment arises 
from the excluded volume interaction. As noted above, 
the force F together with friction from the substrate gives 
rise to velocity matching among the particles. Hence the 
minimal model considered here has the necessary ingre- 
dients to study the collective behavior of self-propelled 
particles. 

Each rod is characterized by the position of its cen- 



ter of mass and its orientation with respect to an exter- 
nal axis, described by a unit vector = (^cos9i,sm9i) 
directed along the axis of the rod. The microscopic dy- 
namics of the z-th rod is controlled by coupled Langevin 
equations, 

(1) 

where Greek indices denote Cartesian components and 
Vij = Ti — Vj . The self-propulsion velocity of magnitude 
Vq is directed along the long axis of the rod. We choose 
coordinates so that the z axis is normal to the plane of the 
substrate. The excluded volume interaction V (vij, 9i, 9j) 
has the Onsagcr form, with 

Vex {^ij,9i,9j) = 1 if particles i, j intersect 

= otherwise. (3) 

This effective interaction can be derived from the colli- 
sion rules that govern the dynamics of 2d hard rods and 
represents a very good approximation to the true mo- 
mentum exchange in the overdamped limit (isj . The dif- 
fusion tensor, D^^p, and the rotational diffusion rate, Dn, 
arc determined by the energy scale associated with the 
excluded volume interaction and the translational and 
rotational damping provided by the substrate. Finally, 
rjia and 77/* are white noise sources with correlations 

{mc.{t)vjpit')) = D^AjHt-t') , (4) 

(^f W^f W) = DnS.,S{t^t') . (5) 

For systems in thermal equilibrium the fluctuation- 
dissipation theorem requires that the amplitudes of the 
correlation of the thermal noise be identical to the diffu- 
sive parameters controlling the excluded volume interac- 
tion, i.e., Dl^p = D\^p and — Dr. This identiflcation 
does not hold for non-equilibrium systems such as the 
one considered here, where the noise has non-thermal 
contributions. For self-propelled suspensions noise arises 
in general from thermal Brownian motion, hydrodynamic 
interactions among the self-propelled units, and intrinsic 
fluctuations in the activity of each self-propelled unit. 
Brownian noise is expected to be very small for particles 
such as bacteria of size larger than a few microns. The 
other two sources of noise are intrinsically non-thermal 
and will generally not obey a fluctuation-dissipation the- 
orem. It can be shown, however, that the equality still 
holds for a dilute solution in the regime of small self 
propulsion velocities [l^. When the self propulsion be- 
comes large the relationship between these diffusion con- 
stants is complicated. In the subsequent analysis, the 
emergent physics of the system is independent of these 
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details. Therefore for simplicity we assume the noise am- 
plitudes are simply given by their equilibrium counter- 
parts, with 



(6) 



For long thin rods the longitudinal and transverse diffu- 
sion constants D|| and D±^ are simply related as D\]^ ~ 
2D± = 2Do, with Dq = kBT/( for a single rod at low 
density where T is the temperature and is a material 
dependent friction constant characterizing the interac- 
tion of the rod and substrate. A similar approximation 
is used for the rotational noise amplitude, with 



QDq 



(7) 



where the last equality holds for infinitely thin rods at 
low density. 

The noise averaged statistical mechanics of a system 
of overdamped Langevin equations is given by a Smolu- 
chowski equation for the A'^-particle phase space density. 
In the low density limit, when two particle correlations 
can be neglected, this results in an efi^ective mean field 
Smoluchowski equation of the familiar form for the one- 
particle probability density, c(r, 0, i), representing the 
density of rods with center of mass at r and orientation 
9 at time t. Its dynamics is governed by the equation 



dc{Y,e,t) 



(8) 



with translational and rotational currents 

''a 

= ~dr ^' r:"' - . (10) 



dV(r,e) ^ dc 

-Dap C \-VoU0C-Dap^ , (9) 



dV{r,0) dc 



In Eqs. ^ and (fTO]) the two-particle Onsager interaction 
defined in Eq. ^ enters in its mean field form 

V(ri,0i) = j dv2 J de2Ve.iru,Oi,02)c{r2,e2,t) 

= y"d^y"d02|Ui XU2|c(ri +1,02,0 (11) 

where ^ = r2 — ri ~ siUi — S2U2 and —£/2 < Si < £/2 
parametrizes the position along the «-th rod measured 
from its center of mass. As shown in Fig. [TJ this is a 
measure of the area excluded by a rod at a point ri, 
oriented in the direction 61 to all other rods in the system. 

We emphasize that the self propulsion mechanism en- 
ters the Langevin equations ([T]) and ^ only in the form 
of a center of mass force. Therefore it gives rise to no mi- 
croscopic torques. This is reflected in the Smoluchowski 
equation above in the fact that the rotational fluxes in 
Eq. (|10p contain no information about the self propulsion. 
The velocity vq enters only the translational flux as an 
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FIG. 1; Geometry of overlap between two rods of length I. 
Here ^ = r2 — ri is the separation of the centers of mass of 
the two rods. 



additional mass flux along the orientation of the rod. In 
the absence of self propulsion (uq = 0), the microscopic 
equations and hence the Smoluchowski equation are both 
invariant when Ui — > — U; for all rods and the system can- 
not exhibit a macroscopic polar state. When ■^o 7^ 0, this 
symmetry is broken at the level of the Langevin equa- 
tions. One might then expect that the system can order 
in a uniform polar state. But since bulk states are con- 
trolled entirely by the rotational fluxes which still possess 
the nematic symmetry, such a homogeneous polar state 
cannot exist. This is not the consequence of any of the 
approximations made to obtain the Smoluchowski equa- 
tion, but is a true property of the non-equilibrium model 
system described by the microscopic dynamics given in 
Eqs. (P) and The self propulsion term does, how- 
ever, affect the fiuctuations the systems, as described in 
the next section. 



III. HYDRODYNAMIC EQUATIONS 

In this section we obtain coarse grained equations ap- 
propriate for describing the dynamics of the system on 
length scales large compared to the length of the rods 
and time scales long compared to the microscopic diffu- 
sion times in the Langevin equations. In this limit we 
expect that the dynamics will be controlled by the con- 
served densities and the fields associated with possible 
broken symmetries. For an overdamped system of the 
type considered here the only conserved quantity is the 
concentration of rods. It can be defined as the zeroth 
moment of the probability distribution c (r, u, t), 



J duc{r,u,t) . (12) 



Both polar and nematic order are in principle possible in 
a gas of self-propelled rods. These can be described by 
a polarization vector P(r, i) and the nematic alignment 
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tensor Qa[3(j, t) defined as first and second moments, re- 
spectively, of c (r, u, t), 



p (r, P (r, t) = J duuc{r,u,t) 



(13) 



p(r,t) Qaf3 (r,t) 



J du (^aUl3 - ^'5a/3^ ^ (r, U, t) . 



^ (14) 

Since each rod has a self propulsion velocity vqU, the 
polarization is also proportional to the self propulsion 
flow field. 

The equations for these continuum fields are obtained 
by taking the corresponding moments of the Smolu- 
chowski equation ([S]) and have the form 



dtp 



V J, 



_d_ 



J, 



Rcx(5 



(15) 



(16) 



(17) 



where 









J ay 






\ Japy / 







{uaUp - \5al3) 



, (18) 



and 



Ra 

Raf) 



du 



{UaUf} - \5al}) ) do ' 



(19) 



with the fluxes J-^ and as given in Eqs. ([9]) and pO|l . 
In order to obtain a closed set of hydrodynamic equa- 
tions, we need to evaluate the constitutive equations for 
the fluxes on the right hand side of Eqs. (|15m7[) . This 
requires various approximations. First, since we are in- 
terested in a long-wavelength description of the system, 
the nonlocal dependence on the concentration field in 
Eq. (fTTj) is expanded in gradients as 

c(ri U2,i) = c(ri,U2) -t- ^q^^i^c (ri, U2) 

+ \^c^p^r,^^r,,c{v^,n2) + 0{V^) ,(20) 



where the expansion has been truncated at second order 
in the gradients. This amounts to a hydrodynamic de- 
scription up through Navier-Stokes order in the coarse 
grained densities. 

Secondly, the balance equations for the polarization 
and the nematic order parameter couple to higher mo- 
ments of the probability distribution. To obtain a closed 
set of equations we need an additional assumption that 
allows these higher moments to be expressed in terms 
of the first three moments of the probability distribu- 
tion. We assume that at long times (times longer than 
all microscopic diffusion times) the higher moments be- 
come functionals of the first three moments [l^ and the 
probability distribution has the form 



c(r,u,t) ^ c(u, [va (r,t)]) 



(21) 



where ya{Y,t) = {p (r, t) , P (r, <) , Qq/j (r, t)} and the 
square brackets denote the fact that the probability dis- 
tribution c is a functional of these fields. Once such a 
functional assumption is made, the Smoluchowski equa- 
tion can be self consistently solved to obtain the prob- 
ability distribution in terms of the slow variables in the 
system and their gradients. In the context of equilib- 
rium systems, this procedure is implemented using slow 
variables that are hydrodynamic, i.e., they are the con- 
served quantities of the system. But as has been shown 
in the context of granular fiuids, it can carried out using 
overdamped variables as in the case at hand as well [2^ . 

Here this procedure is implemented in a much simpler 
context. It is assumed that the probability distribution 
is a linear functional of the slow variables and the solu- 
tion to the Smoluchowski equation is obtained to lowest 
order in the gradients in the slow variables of interest 
here. Then, it follows immediately that the probability 
distribution is given by 



c(r,u,i) = -^p(r,<) l + 2P(r,<)-u 

ZTT 



+4Qq/3 ( UaUp - -;^5af3 



(22) 



Using the two approximations given by Eq. (|20p and 
Eq. p2|) . a set of closed macroscopic equations for the 
density, polarization and nematic order parameter are 
obtained. Including only diffusion and self propulsion 
contributions to the fluxes, these equations are given by 



dtp + V^V ■ (pP) = ^V^p + ^^dadfiipQ^p) 



(23) 



dtipPc) + l^odo^P + vodpipQc^p) = -DrpP^ + ^v\pP^) + ^d^V ■ (pP) (24) 
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dcipPp) + dp{pP^) - So^pV ■ (pP) 



^0 



9a 9/3 



Do 



-~4:DB,pQa0 
d-f{daPQl3f + dppQa-f - SapdaPQai) 
I 



(25) 



The excluded volume contributions to the fluxes are given 
in Appendix A. In Eqs. (|23ll25p . self propulsion generates 
convective-type terms that couple the density and the 
nematic order parameter to the polarization, reflecting 
the fact that self propulsion yields a mass flux in the 
direction of polarization. In the following sections, these 
equations are used to study the possible steady state of 
the systems and their stability. 



A. Homogenoues states 

The bulk states of the system are determined by the 
solutions of the homogeneous hydrodynamic equations. 
Dropping all gradients terms, but including excluded vol- 
ume effects, these are given by 



dtp = 0, (26) 



dtpPc = -DrpP^ + -Dnmp'PpQ^p, (27) 



dtpQaP = -4Di^p(l - 7]op/3)Qa/3 , (28) 

where rjo = 2£'^/7t is the excluded volume of a rod. As 
anticipated, the self propulsion velocity vq does not en- 
ter the homogeneous equations. This is an artifact of 
our modeling of self propulsion solely as a center of mass 
force and hence an effective center of mass velocity in the 
overdamped limit in the microscopic model considered 
here. The homogeneous equations are therefore identical 
to those of an equilibrium collection of overdamped hard 
rods. At low density the system forms an isotropic liquid 
state, with p ^ po, P ^ and Qa/^ = 0. The isotropic 
state becomes unstable for po > Pn = ■^Vo = 37r/2^^ 
where the coefficient of Qap on the right hand side of 
Eq. ([^ changes sign. For po > pN a collection of hard 
rods forms a nematic liquid crystal, characterized by bro- 
ken orientational symmetry along the direction of a unit 
vector n known as the director. The alignment tensor 
takes the form 



iafj S [ UaTip - -Sap 



(29) 



This is the well-known result due to Onsager. In mean- 
field theory the isotropic-nematic transition is continu- 
ous in two dimensions [2lj. To obtain the value S of 
the order parameter one must retain terms cubic in Qap 
in Eq. (|28p which arise from higher order correlations 
neglected here. These yield a term ^ wp'^Q^sQisQap, 



where in the present context w is a phcnomcnologi- 
cal parameter, independent of density. One then finds 

S = -^\J ^ {po/pn — 1) in the nematic state. Finally, 
there is no homogeneous polarized state. 



B. Fluctuations in the isotropic State 

In this section we examine the dynamics of fluctuations 
in the isotropic state for po < PJV to study its stability. 
The only hydrodynamic variable in the isotropic state is 
the concentration of filaments, p. Fluctuations in both 
the polarization and the alignment tensor decay on mi- 
croscopic time scales. However, while fluctuations in the 
alignment tensor couple to the density to higher order 
in the gradients and can be safely neglected, polariza- 
tion fluctuations can qualitatively change the nature of 
the hydrodynamic modes at intermediate length scales. 
We therefore consider the coupled dynamics of fluctua- 
tions of the density Sp{r, t) = p{r, t) — po ^^nd polarization 
SP{r,t) = P(r,t) from their homogeneous values p = po 
and P = 0. It is convenient to expand the fluctuations 
in Fourier components 



^(k,t) 



dre'^ 'Sy (r, t) . 



(30) 



and to write the polarization in terms of its compo- 
nents longitudinal and transverse to k, as P(k, t) = 
k(5P||(k,i)+kj_(5PL(k,t), where k = k/|k| andk_L = zxk 
are unit vectors longitudinal and perpendicular to k. The 
transverse component of the polarization decouples from 
the density and will be neglected below. The linearized 
coupled equations for fluctuations in the density and lon- 
gitudinal polarization, which corresponds to splay defor- 
mations of the polarization field, are given by 



dtSp — —Dk^6p + ikvoPo5P\\ 



(31) 



6p 



dtSPn =-[Dr + Dpk^]dP« + ikvo-^ (32) 

2po 

where in the low density approximation considered here 



D = 



1 



Po 
Pin 



7. 



Dp = ^Do 



(33) 
(34) 



Note that the splay diffusion constant Dp can naturally 
be written as Dp = Ki/Q, with Ki the splay elastic con- 



stant and C, = 
we obtain Ki 



Do 



a friction coefficient. At low density 
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li vq = 0, Eqs. ([3T|) and ((32|) decouple: density fluc- 
tuations decay via diffusion and polarization fluctuations 
are overdamped with decay rate Dr. Conversely, at fi- 
nite vq, for t >> D'^^ and long wavelengths, we can 
assume that polarization fluctuations relax quickly to a 
value determined by the inhomogeneous density field, 
i5P|| ~ ikvQSp/{2Dfipo), which yields traveling density 
waves of speed vq. 

In general we look for wavelike solutions with time de- 
pendence given as linear combinations of terms of the 
form ^ g-^^Cfe)*^ where Ai^(fc) are the dispersion relations 
of the hydrodynamic modes of the system. The modes 
controlling the decay of density and splay fluctuations in 
the isotropic state are given by 



Ai, 



DR + {Dp + D)k^ 



±T^^[DR + {D.p-D)k^Y -2vlk\ (35) 



The sign of the real part of the eigenvalues \v{k) controls 
the linear stability of the homogeneous state. It is easy to 
see that in the isotropic state Re[Xu{k)] < for i/ = 1,2 
and all values of parameters, indicating that the isotropic 
state is linearly stable for all po < Pn- The nature of the 
modes changes, however, from diffusive to propagating 
when the argument of the square root on the right hand 
side of Eq. changes sign. At low density Dp > D 
and the system exhibits propagating density waves for 
■^0 ^ vo{k), with 



voik) 



Dr + {Dp ~ D)k-' 
^/2fc 



(36) 



The crossover from diffusive to propagating density ffuc- 
tuations is displayed qualitatively in Fig. [2l There is a 
lower value Wco = \/2Dr {Dp — D) ^ \/DrDo/2 of vq 
below which the behavior is always diffusive For a fixed 
value of vo > Wco propagating waves exist in a range of 
wavevector fed < fc < kc2, with 



D„-D 



It 



7)2 



^ ■ (37) 



This range widens as vq increases, with lim^jQ^oo ^ci ~ 
and lim„„^oo A;c2 = oo. This behavior resembles closely 
the appearance of propagating sound waves at interme- 
diate wavevectors in a compressible fluid that interacts 
frictionally with a substrate [23|. The polarization here 
plays the same role as the flow velocity in the Navier- 
Stokes equations of Re. [11] and the self-propulsion ef- 
fectively lowers the damping of polarization fluctuations, 
yielding propagating waves. 

Using Do ^ PDr, we obtain Vco ~ Do/i ~ (-Dr. It is 
useful to introduce two characteristic length scales as 




(38) 



propagating 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 



FIG. 2: A qualitative description of the crossover from diffu- 
sive to propagating density fluctuations. The solid line (red 
online) is the boundary for vo{k) given in Eq. p6|l . For a 
fixed value vo > v^o the modes are propagating in a range of 
wavevector kd < k < kc2- 



where (diff and isp are the distances travelled by the 
center of mass of a rod due to diffusion and self propul- 
sion, respectively, during the characteristic time scale for 
rotational diffusion, D^^ . The condition vo » Vco of 
large self propulsion corresponds to £sp » £diff- In this 
limit 



(39) 
(40) 



The modes are propagating for all length scales larger 



than £lff/£. 



« 



'diff 



and smaller than £sp. For any 



flxed value of vo, the modes always become diffusive when 
A; ^ 0. Finally, at high density we expect D > Dp. 
In this case propagating waves will exist in a range of 
wavevector for all nonzero values of vo- Furthermore in- 
creasing the density will result to a strong suppression of 
rotational diffusion from entanglement. In this case prop- 
agating waves will exist for essentially all length scales. 



C. Fluctuations in the Nematic State 

Here we examine the linear stability of the nematic 
state, characterized by a uniform density po > pn, align- 
ment tensor = S{noanoi3 — ^Sap), and P = 0. The 
ordered state is symmetric for iio — Sq. The hydrody- 
namic variables in this case are the density and the direc- 
tor. As for the isotropic state, fluctuations in the polar- 
ization, although overdamped, can qualitatively change 
the system dynamics and ultimately render the uniform 
nematic state unstable. They will therefore be incorpo- 
rated in the analysis below. We assume the system to be 
deep in the nematic state and neglect fluctuations in the 
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magnitude S of the order parameter. We consider small 
fluctuations about the ordered state by letting 

p{r,t) = po + Sp{r,t) , 

P(r,t) = ,5P(r,t) , 

Qai3 (r, t) = Qlp + S [rioaSnj_p (r, t) 

+6n±a (r,t)rio/3] , (41) 

where Sn±^ denotes fluctuations perpendicular to Sq. To 
linear order, this is the only fluctuation in the director 



field that preserves |n| = 1. It is convenient to choose 
a coordinate system with the x axis along Sq, so that 
6n± = 6nyy. 



The fluctuations SP in the polarization are over- 
damped at the rate Dn. For t >> D~^^ we neglect 9t(5P 
and eliminate polarization fluctuations in favor of the 
density and director fields. The resulting hydrodynamic 
equations arc given by 



dtSp = [D.,dl + Di_^\yp + D^Spo (1 + 2a) d,V i_ ■ Sn^ , 
dSn^ ^ (^K^dl + KiVa_Vi_ ■ )(5n 

I — 



§(3 + 2a)a.V,^ 

Hb po 



(42) 
(43) 



with 



and a = 





Do 
4 


[3 + 5 + 2a (1 + S-)] , 


(44) 


Di_ 


_ Do 
4 


[3 - 5* + 2a (1 - 5)] , 


(45) 


Ki 




= ^(3 + a), 


(46) 


D 


- 4- 


a dimensionless parameter. 


As 



expected on the basis of symmetry, only the square of 
the self propulsion velocity enters the equations for the 
nematic. Self propulsion enhances diffusion along the 
direction a; of alignment. It also stiffens both the bend 
and splay elastic constants, K2, and Ki. To analyze the 
stability of the nematic state we expand the fluctuations 
in Fourier modes at wavevector k. Denoting by the 
angle that k makes with Ho, the equations for the Fourier 
amplitudes of the fluctuations are 



dtSp = -Dpp (a, (f)) k^Sp - SpoDpn{a, 4>)k'^5ny 
dtSuy = -Dnpia,(j>)k'^- Dnn{a)k'^5ny , 

I 



(47) 
(48) 



with 



-(3 + a) 



Dpp = [D.:^sm^ (f) + D^cos^ 

Do, 
4 

Do{l + 2a) cos (f) sin 4 

Dnp = — - (3 + 2q;) cos (f> sin < 
8 



(49) 
(50) 
(51) 
(52) 



When (f) = 0, i.e., for wavevectors k parallel to rio, density 
and director fluctuations decouple and the modes are al- 
ways stable. To linear order, this geometry corresponds 
to pure bend fluctuations of the director. The modes 
also decouple for (/> = 7r/2, corresponding to k normal to 
the direction of broken symmetry. To analyze the lin- 
ear stability of the nematic state for arbitrary angles (f> 
we must examine the the dispersion relations of the hy- 



drodynamic modes associated with Eqs. (|47p and 
These are easily obtained as 



-2 r 



-(D 



pp 



D„ 



± \/{Dpp - Dnn)^ + ADpnDnp 



(53) 

The eigenvalues are always real, corresponding to diffu- 
sive modes. The mode A_ is always negative, indicating 
stable decay of fluctuations. The mode A+ changes sign 
for 



(54) 



It is easy to verify that this condition can never be sat- 
isfied when vq = 0, i.e., the nematic state is stable in the 
absence of self propulsion. On the other hand, the ho- 
mogeneous nematic state becomes unstable for all values 
of k at large enough dq. The condition for the onset of 
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FIG. 3: (color online) A schematic plot of the boundary 
Vc{(l>,S) above which the homogeneous nematic state is lin- 
early unstable as a function of cp for S — 1 (solid red line) and 
S = 0.75 (dashed blue line). For a fixed value of vq > Vc{S) 
the nematic state is unstable for (pci 4'c2- 



the instability obtained from Eq. yields a boundary 
Vc^fj): S). The nematic state is unstable for vq > Vc{4>i S), 
as shown in Fig. ([3]). For a fixed value of S, there is a 
minimum threshold value Vc{S) required to destabilize 
the nematic state. For a fixed value of vq > Vc{S) the ne- 
matic state is unstable for all spatial gradients such that 
0ci < < 4'c2- For S* = 1, we find limug^oo 0ci = j and 
Imiyg^oo (f)c2 ~ and the system becomes unstable to 
all spatial gradients such that the splay wins over bend. 
Finally, note that the unstable region shrinks in size as 
the magnitude of nematic order in the system decreases. 

The instability of the nematic state arises from a subtle 
interplay of splay and bend deformations and diffusion 
longitudinal and transverse to the direction of nematic 
order. It only occurs for angles (f) > 7r/4, corresponding 
to situations where splay deformations of the nematic di- 
rector exceed bend deformations. It does not, however, 
occur for pure splay (</> = 7r/2), other than at uq — > oo. 
In a nematic, curvature inhomogeneities of the director 
yield mass currents, with splay deformations ~ dySny 
yielding current along the direction x of mean order and 
bend deformations ^ dxSuy yielding current along the y 
direction. Self propulsion enhances both contributions. 
In addition, it enhances mass diffusion both longitudinal 
and transverse to the direction of order, as indicated in 
Eqs. ([13]) and However, while the enhancement of 

Dx from self propulsion grows as the nematic order in- 
creases and is maximum at S" = 1, the enhancement of 
D± drops as one moves deeper into the nematic state 
and vanishes when S = 1. Thus, if a curvature in the 
director field induces a mass flux that requires relaxation 
through appreciable diffusion perpendicular to the mean 
director, then, beyond the threshold value Vc, the sys- 
tem becomes unstable. A fluctuation in the director field 
corresponding to a large splay and a small bend results 
in exactly such mass fluxes. Conversely, when bend ex- 
ceeds splay longitudinal diffusion is the main mechanism 
that restores orientational order. This is enhanced by 
self-propulsion and the system remains linearly stable. 



The analysis described above only includes terms up 
to quadratic order in the gradients and predicts that the 
nematic state is unstable on all length scales. We expect 
that terms of higher order in the gradients arising from 
excluded volume and diffusion currents will stabilize the 
nematic state at large wavevectors [l^], setting a length 
scale for the typical size of the ordered regions, not unlike 
what has been observed in the simulations by Peruani et 
al. 
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IV. SUMMARY AND DISCUSSION 

In this paper we have discussed the collective dynamics 
of self propelled hard rods on a frictional substrate. The 
rods interact through excluded volume interactions. The 
self propulsion was implemented as a center of mass force 
acting on each rod that propels it along its long axis. Self- 
propulsion breaks the nematic symmetry at the level of 
the microscopic equations of motion, but since it yields 
only a center of mass force on each rod it is insufficient to 
generate a macroscopic polarized state. The existence of 
a polar state requires either microscopic torques that turn 
the particle towards an externally determined direction 
(for instance as seen in the case of chemotaxis in bacterial 
motion) or an aligning interaction that can distinguish 
the two ends of the extended object, which a physical 
excluded volume interaction of rods does not. As a result, 
the only possible bulk states of the self-propelled systems 
are those possible in equilibrium : isotropic and nematic. 

Self-propulsion, does, however have a profound effect 
on the nature of the fluctuations in each of this states. 
In the isotropic state self-propulsion yields the appear- 
ance of propagating waves in a range of wavevectors, as 
shown in Fig. [21 This result is closely analogous to the 
appearance of sound waves in a compressible fluid on 
a frictional substrate, predicted many years ago by Ra- 
maswamy and Mazenko [23 |. This phenomenon has not 
been observed in equilibrium systems where it requires 
a very small value of the friction. In self-propelled sys- 
tems, self-propulsion itself effectively lowers the friction 
and should yield a wide range of parameters where prop- 
agating density waves may be observable. Propagating 
density waves may indeed have been seen in the collec- 
tive dynamics of epithelial cells on a substrate [1^ . In the 
nematic state fluctuations in the local polarization yield 
mass fluxes and can destabilize the homogeneous nematic 
state. This has been previously observed in the numerical 
study of [ij]. Further, before the onset of the instability, 
these same director fluctuations and the resulting anoma- 
lous mass flux associated with them account for the large 
number fluctuations predicted in [l0| . 

Our work has several limitations. First, we neglect all 
correlations in deriving the Smoluchowski equation and 
in the functional assumption made to obtain the hydro- 
dynamic equations. On the other hand, our hydrody- 
namic equations have precisely the structure predicted 
on the basis of pure symmetry considerations. This sug- 
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gests that the behavior obtained here may be generic. 
Therefore, further insight into the transport processes 
can be obtained from numerical simulation of the model 
system and comparing the results to the inherently low 
density theory presented here. Secondly, the overdamped 
Langevin microdynamics that describes the hard core in- 
teraction as an equilibrium mean-field excluded volume 
effect may not be an adequate starting point to incorpo- 
rate self-propulsion. It can be shown that the momen- 
tum transfer that occurs in the interactions between self- 
propelled hard rods modifies the diffusion processes in 
the system [l^. These effects can be captured by start- 
ing with the true Langevin equations describing a system 
of hard rods on a substrate and then obtaining the asso- 
ciated Fokker-Planck equation, as will be described in a 
future publication [i8| . 

V. ACKNOWLEDGEMENT 

Research supported by the NSF through grants DMR- 
0305407 and DMR-0705105. MCM also received support 

I 



re: 



jex 



jex 



where the notation 

^a/Sfa — SapSja + Sa^ySflcr + SaaSpj (A4) 
I 



from the Institut Curie in Paris via a Rotschild-Yvette- 
Mayent sabbatical fellowship. She thanks both the In- 
stitut Curie and ESPCI for their hospitality during the 
completion of some of this work. Finally, we thank Sri- 
ram Ramaswamy for useful discussions and for alerting 
us to Ref. m. 



APPENDIX A: HYDRODYNAMIC FLUXES 



In this appendix, we give the expression for the ex- 
cluded volume contributions to the translational and ro- 
tational fluxes defined in Eqs. (fT8|) and (fT9|) . These 
have been obtained assuming a low moments closure for 
the one-particle distribution function of the form given 
in Eq. ([2T|) . The excluded volume contributions to the 
translational fluxes are given by 



(Al) 



(A2) 



(A3) 



has been introduced for compactness. The excluded vol- 
ume contributions to the rotational fluxes are 



1 / 3 \ 2 7 2 

--^DqTIo {Aapjcr + 'iSaaSpj) pP^d^p + —DqTIo {AafSuj-, + GSaujSpj) pPad^pQa-i 
o io 

+ -r£'o??o[(AQ/3aJ7 - SapS^-^) pPuidcrpQa-y + pPujd^pQap] , 

io 

1 3 

lZ o 



1 

3' 

— P {7S^aQaP + SafQtSa + SpaQaj + SpjQaa + SaaQpj — "^Sa^Qjcr) d^P 

+P^ {SaaQl3LU + SpcrQaiu — ^a/lQauj) Q-y^ + P^ {^ajQ/SLU + Sfj^Qauj — ^afjQ-yu) Qaui 
+ ^P^ {SajSp^ + &o,a5f)^ - l^5^p5^„) Tr (Q^) + Qp^S^aQlp] , 

I 



-DRrio — [pP^d^df3 {pQap) + {pPfidad^ - pPadpd^) {pQp^)] , (A5) 
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riex 

-n-a/3 — 



(A6) 
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